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Abstract 

We study the non-monotonic force-extension behaviour of hehcal ribbons using a new model 
for inextensible elastic strips. Unlike previous rod models our model predicts hysteresis behaviour 
for low-pitch ribbons of arbitrary material properties. Associated with it is a first-order transition 
between two different helical states as observed in experiments with cholesterol ribbons. Numerical 
solutions show non-uniform uncoiling with hysteresis also occurring under controlled tension. They 
furthermore reveal a new uncoiling scenario in which a ribbon of very low pitch shears under 
tension and successively releases a sequence of almost planar loops. Our results may be relevant 
for nanoscale devices such as force probes. 

PACS numbers: 46.32.-Kx, 46.25.Cc, 87.10.Pq 
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Helical ribbons are common in nature Q, Q] and hold great promise as parts of nano- 
scale force probe devices j^. Smith et al. jjstudied the force-response behaviour of helical 
ribbons that appear as metastable intermediates in the process of cholesterol crystallisation 
in the native gallbladder bile. Cantilever tests showed linear, Hookean, and reversible be- 
haviour at low tension, while at sufficiently high tension the ribbon becomes metastable and 
separates into a straight and a helical domain. To explain this behaviour they proposed a 
phenomenological model for an inextensible strip. Intrinsic curvature entered the elastic en- 
ergy through an isotropic term that can be interpreted as describing the difference in surface 
tension between the ribbon's inner and outer surfaces. The model is only valid for helical 
configurations and the spontaneous helical pitch of the ribbon is determined by anisotropy 
of the strain tensor, rather than by intrinsic geometry. The model shows multistability with 
large and small pitch solutions, and associated hysteresis behaviour. However, recent X-ray 
diffraction studies in show that these oho.este.0. nbbons are single crystals and suggest 
that the crystalline structure rather than the surfactant layers determine the helical form. 

Kessler & Rabin considered an elastic rod model and by direct energy minimisation 
were able to produce non-reversible (hysteresis) force-extension behaviour for helical springs 
of small pitch and large ratio of torsional to bending stiffness. This behaviour was con- 
firmed analytically by Zhou et al. who used a rod model that also allowed for transverse 
anisotropy. A single hysteresis cycle was found for large stiffness ratios. In [7] dynamical 
simulations were performed using an extensible isotropic rod model revealing the possibility 
of several abrupt changes in the force response of a stretched helix, again in the case of 
unusually large stiffness ratios (corresponding to negative Poisson ratios). 

We propose a new model of an elastic inextensible strip. Such a strip can only bend; 
thus, if its unstressed surface is developable, i.e., has vanishing Gaussian curvature, then it 
will remain so after deformation. We believe that this model is better suited for crystalline 
structures. We will show that it predicts hysteresis behaviour for realistic values of Poisson's 
ratio and without requiring material anisotropy. 

Consider an inextensible ribbon that, when developed onto a plane, makes a strip of 
width 2w bounded by two parallel straight lines. By the property of developability, the strip, 
however deformed, can be reconstructed from its centreline r{s)^ which itself is determined, 
up to Euclidean motion, by its curvature k(s) and torsion t{s). Here s is the arclength along 
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FIG. 1: Left: A developable surface is made up of straight generators at angle (3 to the tangent t to 
the centreline r; b is the binormal. Right: A developable helical ribbon under tensile end loading. 

the centreline. More precisely, we can parametrise the strip as 

x{s,t) = r{s)+t[h{s)+r^{s)t{s)], (1) 
r(s) = ?7(s)k(s), s = [0, L], t= w], 

where t is the unit tangent vector and h the unit binormal. The generator of the developable 
surface of the strip makes an angle /3 = arctan(l/r7) with the positive tangent direction 
(Fig. DP). 

Let N be the unit normal to a smooth surface and S{T) = —OtN the shape operator 
acting on a unit tangent T to the surface. The bending energy for a Kir chhoff- Love shell of 



thickness 2h can be written as the following integral over the surface 
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U = ^ jjwMAS))' + (1 - ^) trUASnda, (2) 

where D = 2E / [3{1 — u"^)] is the fiexural rigidity, u is Poisson's ratio, E is Young's modulus 
and AS* = 5* — 5*0; here and later the subscript refers to the undeformed state. Equation ([2]) 
can be rewritten in the form llC 



U = D J J [{1 + u){AHy + {1- iy){{AAy +4AAo sin' e)]da, (3) 
where AH = H — Hq, H is the mean curvature, 6 is the angle between the principal 



curvature axes in the deformed and undeformed states, and A = H'^ — G, G being the 
Gaussian curvature. 



For developable surfaces, H = A = k,i/2 and K = (ki and K2 = are the principal 
curvatures). We shall assume the undeformed strip to be helical, with constant pitch angle 
7r/2 — Pq, and to lie on the surface of a cylinder of radius Rq (Fig. ^^ight). Then Eq. ([3]) 
becomes 

D 

U = - / [kI- 2Ki/ti,o(l + sin^ 6) + /t? q] dtds, (4) 

^ Jo J-w 

with Ki^Q = I/Rq and 6 = P — Pq. On using Eq. ([21) the t-integration can be carried out 
analytically, as in the case without intrinsic curvature (l2|], and we arrive at 



with 



U = Dw / g{K, T], r]')ds + const (5) 
Jo 



2 2 

g{K,r],r]') = {l + rf) V [wr]') - — k [l + + q{ri - rj^f] , (6) 

/to 



where g = (z/ — 1)/(1 + t^q), 770 = cot /?o and the prime denotes the derivative with respect to 
s. Note that for strips with intrinsic curvature (-Rq < 00), equilibrium shapes will depend 
on the material properties through Poisson's ratio. Also note that in the limit of narrow 
strips, wrj' —>■ 0, we have V{wri') 1 and no derivative enters the integrand in Eq. 

A different model of a thin and narrow elastic strip may be obtained as a limiting case 
of the thin anisotropic elastic rod by pushing one of the principal bending stiffnesses to 



infinity 



131 ]. In this model the material frame of the rod coincides with the Frenet frame of 



its centreline. We shall call it a Frenet rod. Its surface is swept out by the binormal of the 
centreline and is not obliged to deform isometrically. If the centreline has intrinsic curvature 

n 

kq and intrinsic torsion tq then the elastic energy reads [5| 

f/F = f («: - ^0? + jir - Tofds, (7) 

where B and C are the bending and torsional stiffnesses, respectively. Assuming a rectan- 
gular cross-section, we have B = AEwh^/3, B/C = (l + z/)/2 and Eq. (JTj) may be rewritten 

as 

B 

Up = — / gF{K,, 7]) ds + const (8) 

2 Jn 
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with 



For both the strip and the rod model enerery minimisation is a ID variational problem 

n 

cast in a form that is invariant under Euclidean motions. Following [1^ , the Euler-Lagrange 
equations can be immediately written down in the form of six balance equations for the 
components of the internal force F and moment M in the directions of the Frenet frame of 
tangent, principal normal and binormal, F = {Ft, Fn, Ff,)'^ , M = {Mt, Mn, Mb)'^ , and two 
scalar equations: 

F' + u; X F = 0, M' + ujxM + txF = 0, (10) 

d^g + r]Mt + Mfe = 0, {d^,g)' - d,g - nMt = 0, (11) 

where a; = n{r], 0, 1)"^ is the curvature vector. The equations have |-Fp and F ■ M as first 
integrals. We note that in the variables (k, 77, 77') the first equation in Eq. (ITT!) is algebraic 
and the second equation becomes so in the limit of a narrow strip {w = 0) and for the rod 
(Eq. (ED). 

We first consider helical solutions, i.e., we set k, = const 7^ and r = const. In addition, 
we assume that F and M are constant in the local (Frenet) frame. Then Eq. ffTOj) implies 

F„ = 0, M„ = and 

K{Mt - vM,) = Ffe, Ft = vFb. (12) 

Thus, the vector F is directed along the helical axis (z) and Ft = F cos (3, Fh = F sin (3. The 
magnitude of the force can be obtained from Eqs. (fTTj) . ([T2|) : F = -^/l + [nrjdi^g — (1 + 
rf)dr^g] (negative for tension). 

We assume that no moment acts about the helical axis, i.e., the moment is orthogonal 
to the force: F ■ M = 0. Then the first Eq. (fTTl) simplifies to d^g = 0. This allows us to 
eliminate the curvature and to obtain an explicit expression for the force as a function of 
only one variable, for which it is convenient to choose the angle (3. We furthermore introduce 



the normalised extension ( = z{L)/L = cos/? so that rj = (/ ^yl — C^. For our functional g 
we find for the normalised force / = FRq 

f = ^7^[2(1 + ^) sin2(/3 - /?o) + (1 - ^) sin4(/3 - /3o)], (13) 
^ Sin LJ 

while the analogous expression for the Frenet rod functional gp is 

, _ ^[(1 + I/) sin/?osin/? + 2cos/3ocos/?] sin(/3 - /?o) /^^^s 
sin/3[(l + t/)sin2/3 + 2cos2 /3Y ' 
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arctan(2) 



FIG. 2: The I'-Po parameter plane for helical solutions with regions of different force-extension 
curves Eq. ([13]) indicated. The dashed curve is given hy u = — cot^ Pq. Dark shaded is the physical 
region of multistability. Light shaded is the region, bounded by the curve given by cos Po = 



\u\l \J\ + (1 — z^)^, where the helical solutions are stiffer for the strip than for the rod. /?* = 
arctan (3/2) 

This expression coincides with the corresponding formula for the rod with circular cross- 
section in [3] after replacing the torsional-to-bending stiffness ratio in the latter with that 
for a rectangular cross-section. 

For small deformations the response of the helical ribbon may be approximated as 



Hookean (linear) with spring constant (9^/1 



IC=Co 



4(1 — I/) sin ^/5o. The corresponding 



value for the Frenet rod, | def y 



ion 0,[li 



iC=Co 



4[1 + + (1 — z/) cos^/3o]~"^, may be obtained from 
the classical expression [7|, Il5| by replacing the torsional-to-bending stiffness ratio with that 
for a rod of rectangular cross-section. Which of the two is stiffer depends on parameters, as 
indicated by the light shading in Fig. [2j 

It is straightforward to show that for positive v and for arbitrary /^o the force / vanishes 
only at two extensions: C = Co (undeformed state) and C = Ci < 0, Co + Ci = 1 (helical axis 
orthogonal to the undeformed flat direction), while for negative f has 4 roots. Fig. [2] shows 
force-extension curves Eq. f[T^ as insets in the z/-/5o parameter plane. Non-monotonicity of 
these curves (hysteresis) implies the coexistence of two stable solutions at given force /, one 
of high and one of low pitch. A first-order phase transition between specific values of these 
pitches occurs at a certain force for which the total energy densities g + F( of the two helices 
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FIG. 3: Force-extension diagram for the inextensible strip (n = 10, i/ = 0.3, Pq = 85°). The dashed 
curve corresponds to Eq. ()13p for hehcal solutions. 

('phases') are equal. Solutions may then be constructed composed of arbitrary pieces of each 
lelix. Such solutions describe the phase-separated configurations observed experimentally in 
The domain of hysteresis or first-order transitions in the iy-(3o parameter plane is bounded 
by critical curves {u = 0} and 7 of second-order phase transitions or cusps given by i9^(/ = 0. 
Note that the maximum physically realisable value of Poisson's ratio, u = 1/2, ensures an 
interval of initial pitches with phase transitions. The curve in Fig. [2] is the critical curve 
for the Frenet rod [sl, Eq. (fl^ . It bounds a much smaller region of multistability 17]. 

Figs [3H6] show results of numerical calculations in which we solved Eqs ffTOl) and ffTTj) 
for a finite-length and zero-width ribbon of n helical turns subject to zero-moment {M = 
0) boundary conditions (thus emulating the conditions in the experiments in I2I). Fig. [3] 
gives the response curve and superimposed ribbon shapes showing non-uniform uncoiling. 
Note that ( is now the end-to-end distance. For extensions less than about 0.8 the curve 
closely follows the (dashed) helical curve. At larger extensions more complicated hysteresis 
behaviour is found, with jumps occurring not only under controlled force (known from j^), 
but also under controlled extension. This multistabihty allows one again to construct multi- 
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FIG. 4: Force-extension diagram for the inextensible strip: n = 5 (solid), n = 10 (dashed); v = 0.3. 
Left: Po = 85°. Right: Pq = 45°. Dotted curves are for exact hehcal solutions. 
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FIG. 5: Force-extension diagram for the Frenet rod. Left: n = 10, = 0.1, /?o = 88.24°. Right: 
n = 5, = 0.3, Po = 85°. Dashed curves are for exact helical solutions. 

phase solutions. Such solutions will in general not be smooth at the interface and, in order 
to be supported mechanically, would require special joints. Smooth solutions may however 
exist for certain parameters and boundary conditions. 

Fig. m compares results for various values of n and Pq showing that non-monotonic force 
response is only observed for sufficiently large (3q, a result also found for the Frenet rod 
model j^. The integer n controls the number of 'bumps'. As cx3 the helical curve is 
approached. Fig. ^Left also shows, through the slopes of the curves at small (, that a helix 
is stiffer in tension than a real ribbon that is allowed to deform from its helical state, as is to 
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FIG. 6: Force-extension curves for 3-, 4- and 5-turn helices of very low pitch (/?o = 89°). Left 
V = 0.2. Right: u = 0.05. Unlike in previous figures the force here has been normalised by the 
common factor Rq = l/(27rn*)^ with = 5. Dashed curves are for the planar elastica for n = 2 
(resp. 3) to 5. For small C, these curves have d(^f\^^^^ = — |(n/n=K)^. 

be expected given the helical constraint. For comparison, Fig. Ogives results for the Frenet 
rod. At large C, the force response is quite different from that of the inextensible strip, which 
shows larger oscillations. 

Fig. E] shows curves at a pitch angle close to 0° (at which value the centreline of the ribbon 
becomes planar and is described by Euler's elastica with intrinsic curvature). An interesting 
unlooping scenario is observed near this limit with the first response of the helix being a 
shearing to a (near) planar shape. The subsequent force-response is then essentially that of 
the planar multi-looped noninfiectional elastica with zero end moment whose response curves 
are included in dashed lines (the remarkable agreement worsens slightly for larger n since the 
loops are not exactly in one plane). 3D unlooping transitions (again with hysteresis cycles) 
occur in which successive loops are lost one by one. For smaller Poisson ratio the hockles 
get tighter. No self-intersections are encountered in all our computations with w = 0. 

We have explained the tension-induced multistability of helical ribbons by means of a 
geometrically-exact model of an inextensible strip for which we derived new equilibrium 
equations. Two different types of non-monotonic force-extension behaviour have been found, 
one approximated by the unwinding of an exact helix, the other by transitions between curves 
of the looped planar Euler elastica. Unlike in previous rod models we observe this behaviour 
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without the need to choose very large torsional-to-bending stiffness ratios, which are hard 
to justify in the frame of continuum elasticity theory. In our model the torsional stiffness 
arises naturally as a consequence of the geometric constraint of developability of the ribbon's 
surface on account of the crystalline structure of the material. The predictions of this model 
could conceivably be verified in the controlled environment of nanoribbons [3]. Alternative 
boundary conditions could easily be accommodated, as could a finite width w, while the 
model could be extended to allow for material anisotropy. The multiple hysteresis behaviour 
in helical ribbons could be exploited in the design of force probes or bi- or multistable 
nanoswitches 161 ]. 

This work was supported by the UK's Engineering and Physical Sciences Research Council 
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